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ABSTRACT 


An ever increasing concern on the stratosphere ozone depletion has led to a world wide 
ban on chlorofluorocarbons (CFCs) and prompted rigorous search for alternative 
refiigerants with zero Ozone-Depletion Potential (ODP) and lower Global Warming 
Potential (GWP). To develop an optimum alternative to replace an existing CFG pure 
fluid or mixture requires accurate thermodynamic information (especially vapor-liquid 
equilibrium behavior) of not only the old refrigerants, but also of the possible 
replacements. In the present study, the capabilities of the Wong-Sandler (WS) model 
using different G^ equations (NRTL. UNIQUAC, UNIQUAC-R) have been explored in 
the equation of state (EOS) modeling of isothermal binary vapor-liquid equilibria (VLE) 
for a vastly extended range of refrigerant mixtures and conditions. The results showed 
that the WS mixing rule using any of the three G^ equations in all cases could be 
accurately fitted to binary refrigerant VLE data. Similar results obtained with the 
UNIQUAC-R and UNIQUAC models strongly suggests that the entropic part 
(combinatorial term) has little influence on the accuracy of VLE predictions over the 
range of conditions studied. The WS model showed good convergence ability, which is 
of vital importance in VLE calculations, due to its robustness to parameter variation. Its 
accurate extrapolation capabilities would be useful in extending limited experimental 
data, especially for newer and unconventional refrigerant mixtures. 
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CHAPTER 1 


INTRODUCTION 


In accord with an increasing concern on the stratospheric ozone depletion, the timetable 
toward a complete ban on the production and consumption of chlorofluorocarbons 
(CFCs) and hydrofluorocarbons (HCFCs) has been accelerated to phase them out. The 
study of the ‘ozone hole’ above the Antarctic established the relationship between the 
chlorine and ozone content in the stratospheric zone and dictated that the chlorine oxide 
(CIO') radical is mainly responsible for the ozone destruction. This radical is formed by 
chlorine atoms released by chlorofluorocarbons (CFCs) and hydrochlorofluorocarbons 
(HCFCs), used mainly in the refrigeration industry. The environmental regulations that 
have been imposed on all countries in order to tackle the problem of ozone depletion, 
forced the refrigeration industry to totally phase out CFCs in 1995, and HCFCs in the 
following years (2015 for Western Europe) (Clodic and Sauer, 1994). The recommended 
worldwide ban on many low-molar-mass chlorofluorocarbons has prompted rigorous 
search for alternatives with zero Ozone-Depletion Potential (ODP) and lower Global 
Warming Potential (GWP). 

The chlorofluorocarbons were primarily used as refrigerants, but are also imponant as 
blowing agents in polymer foam manufacture and cleaning solvents in electronic circuit 
manufacture. In many cases, such as refrigeration applications (which is most important 
area of use of CFCs), it is necessary to develop alternative fluids with thermodmamic 
properties which are similar to CFCs they are to replace so that existing equipment can be 
used with a minimum of modification. For this purpose first hydro fluorocarbons (MFCs) 
and hydrochlorofluorocarbons (HCFCs) and their mixmres were suggested, though now 
fluorinated alcohols and their mixtures with HFCs and HCFCs and other fluids are also 
being considered (Sauermann et ah, 1993; Laugier et ah, 1994; Nowaczyk and Steimle, 
1992). The use of mixtures can result in energetic improvements, extensions of the 
application range and favourable influence to undesirable properties of one of the 
components, such as flammability, insolubility in oil or green house effect implications. 



For example, whereas HFC-134a is one of the promising candidates as an alternative 
refrigerant because it contains the hydrogen atom and has short atmospheric live resulting 
in low ODP and GWP, there are disadvantages in using pure HFC-134a. It has poor 
compatibility with lubricants due to the absence of chlorine atoms. This disadvantage can 
be eliminated by mixtures with other refrigerants (e.g. HCFC-124) which have a good 
compatibility with lubricants. Azeotropic mixtures have merit when used as refrigerants 
because they show behavior similar to pure components. Thus the list of possible choices 
for new refrigerant mixtures to replace the CFCs is very large and growing. To develop 
an optimum alternative to replace an existing CFC pure fluid or mixmre requires accuraie 
thermodynamic information (especially vapor-liquid equilibrium behavior) of not only 
the old refrigerants, but also of the possible replacements. To obtain detailed 
experimental data for all promising substitutes is expensive and time consuming. Ideally, 
a predictive model such as UXIFAC would be useful in such cases, at least to screen 
possible candidate mixtures for their suitability in a specific application. How'ever, recent 
attempts (Klieber, 1994) indicate that current UNIFAC group definitions and interaction 
parameters for such fluids are either non-existent or do not lead to accurate predictions. 

At present, the best choice for the design engineer is to collect a minimum amount of 
data at select temperatures and pressures, and then to use these data with a model to 
extrapolate phase behavior to other temperatures and pressures. The vapor-liquid 
equilibrium modeling of refrigerant mixtures can be done using either activity coefficient 
(the so-called gamma-phi) models or equations of state (EOS). Usually EOS modeling is 
preferable because in refrigeration design not only is VLE important, but so are other 
thermodynamic properties (such as density and enthalpy) and these can also be obtained 
from the same EOS model. 

In order to accurately correlate and'or predict VLE with equations of state, it is necessary' 
to have an EOS that can represent pure component saturation pressures accurately. There 
have been significant improvements in the accuracy of cubic EOS for the prediction of 
pure component properties in the last two decades (Sandler et al., 1994). Much effon in 
recent years has been devoted toward developing multiparameter mixing rules. 



especially those that incorporate the excess free energy ( or A^) to extend the cubic 
EOS for the phase equilibrium calculations of complex mixtures (Orbey and Sandler, 
1995a; Heidemann, 1996). There is a rapidly growing literature exploring techniques for 
implanting specific behavior into the cubic equation of state through the mixing rales. 
Models that combine equations of state and G^ models can be categorized into two 
groups: those that make this link at infinite pressure (Huron and \Tdal, 1979; Wong and 
Sandler, 1992), and those that make this link at low or zero pressure (Ivlollemp, 1986; 
Heidemann and Kokal, 1990; Michelsen, 1990a.b, and Dahl and Michelsen, 1990). 

Attempts to incorporate excess Gibbs energy model information into equations of state 
were initiated by Huron and Vidal (1979) who matched the G^ of the solution model to 
that of the EOS at the state of infinite pressure. However, their mixing rule has not 
become widely used because the available G^ parameters at low pressures can not be 
used in their mixing rules. Wong and Sandler (1992) and Wong et al. (1992) incorporated 
the second virial coefficient as well as solution model into the mixing rule obtained by 
identifying the excess Helmholtz free energy (A^) of a solution model into that of an EOS 
at infinite pressure. Available low-pressure solution model parameters are found :o be 
useful in this mixing rule because of the insensitivity of A^ to pressure. This theorerically 
correct mixing rule provides a unified and consistent way of combining two-parameter 
cubic equations of state and excess free energy models to describe the phase behaMor of 
a wide variety of mixtures. The Wong-Sandler model has a broad range of applicability; 
it can be used as a multiparameter correlative model to accurately fit phase equilibrium 
data of highly non ideal mixtures (Wong and Sandler, 1992), or it can be used as a 
predictive model (Wong et ah, 1992, Orbey et ah, 1993) with parameters obtained from 
G^ information available in data banks. 

In this work, the capabilities of the Wong-Sandler (WS) model have been explored in the 
EOS modeling of isothermal binary vapor-liquid equilibria for a vastly extended range of 
refiigerant mixtures and conditions. Following Deiters (1997), Table 1 presenis .the 
engineering nomenclature (‘R nomenclature’) as defined by ASECR-\E (American Society 
of Heating, Refrigeration, and Air Conditioning Engineers), the chemical formula, and 



Table 1: Refrigerant nomenclature 


S.N. 

R Nomenclature 

Chemical formula 

lUPAC name 

1 

R-12 

CFCI 2 H 

Dichlorofluoromethane 

2 

R-13 

CCIF 3 

Chlorotrifluoromethane 

3 

R-22 

CF 2 CIH 

Chlorodifluoromethane 

4 

R-23 

CF 3 H 

Trifluoromethane 

5 

R-32 

CF 2 H 2 

Difiuoromethane 

6 

R-40 

CCIH 3 

Chloromethane 

7 

R-13B1 

CF 3 Br 

Monobromotrifluoromethane 

8 

R-114 

CF 2 CICF 2 CI 

1 ,2-Dichloro-1 ,1 ,2,2-Tetrafluoroethane 

9 

R-116 

CF 3 CF 3 

Perfluoroethane 

10 

R-123 

CFCI 2 CF 2 H 

Dichlorotrifluoroethane 

11 

R-124 

CFCIHCF 3 

1 -Chloro-1 ,2,2,2-T etrafluoroethane 

12 

R-125 

CHF 2 CF 3 

pentafluoroethane 

13 

R-134a 

CF 3 CFH 2 

1 , 1 , 1 , 2 -Tetrafluoroethane 

14 

R-140a 

CCI 3 CH 3 

1 , 1 , 1 -T richloroethane 

15 

R-141b 

CFCI 2 CH 3 

1 ,1-Dichloro-1-Fluoroethane 

16 

R-142b 

CF 2 CICH 3 

1 -Chioro-1 , 1 -Difluoroethane 

17 

R-143a 

CF 3 CH 3 

1 , 1 , 1 -T rifluoroethane 

18 

R-150a 

CCI 2 HCH 3 

1,1-Dichloroethane 

19 

R-152a 

CF 2 HCH 3 

1,1 -Difluoroethane 

20 

R-236fa 

CF 3 CF 3 

1 ,1 ,1 ,3,3,3-Hexafiuoropropane 

21 

R-290 

CH 3 CH 2 CH 3 

Propane 

22 

R-600a 

C 4 H 10 

Isobutane 
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the lUPAC name of all the refrigerants involved in the 37 systems of binary refrigerant 
mixtures of CFCs, HCFCs, MFCs, and HCs studied. Liquid phase (activity 
coefficient) models incorporated in the EOS/G^ approach are the commonly used NRTL, 
UNIQUAC and UNTQUAC-R models, where the UNIQUAC-R refers to the activity 
coefficient model that incorporates only the residual part of the UNIQUAC model. It is 
recognized that, attractive forces are related only to the residual part of the UNIQUAC 
model (Brandani et al., 1998). The results indicate that the combinatorial part of the 
UNIQUAC model has little influence on the accuracy of phase equilibrium calculations 
in the pressure range of interest. 



CHAPTER 2 


LITERATURE REVIEW 

Owing to their ozone-depleting properties, chlorofluorocarbons (CFCs) are now regulated 
under the Montreal Protocol (Kruse and Hesse, 1988). Several hydrofluorocarbons 
(HFCs) are now considered as alternative refrigerants that are harmless or less harmful to 
the ozone layer. These refngerants usually form a nonazeotropic mixture, which has 
recently been considered as the prospective working fluid for refrigeration and heat pump 
systems resulting in improved coefficient of performance (Kruse, 1981). In order to select 
the optimal mixture composition for design and operation, it is necessary to know the 
basic thermodynamic data, especially vapor-liquid equilibria. However, only a few 
systems have been investigated up to now. At present, the best choice for the design 
engineer is to collect a minimum amount of \'T.E data at selected temperatures and 
pressures and then use these data with a model to extrapolate phase behavior to other 
temperatures and pressures. 

The vapor-liquid equilibrium modeling of refiigerant mixtures can be done using either 
activity coefficient (the so-called gamma-phi) models or equations of state (EOS). 
Usually EOS modeling is preferable because in refiigeration design not only is XTE 
important, but so are other thermodynamic propenies (such as density and enthalpy) and 
these can be obtained from the same EOS model. In order to accurately correlate and or 

predict VLE with equations of state, it is necessary' to have an EOS that can represent 

\ 

pure component saturation pressures correctly. In recent years several such equations 
have been suggested (Sandler et ah, 1994). Cubic equations of state (EOS) are widely- 
used for thermodynamic property modeling in industrial design, and have proven to be 
extremely effective (Mathias and Klotz, 1994). They allow a consistent description of 
phase behavior for mixtures. The extension of a cubic, two-parameter EOS to mixtures 
requires mixing rules for the parameters. For relati\'ely simple mixtures, it is common to 
use the van der Waals (vdW) one-fluid mixing rules. VLE data for the CFCs are 
relatively easy to correlate with cubic equations of state using the conventional one- 
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parameter van der Waals (vdW) mixing rules. Also, at least in the narrow ranges of 
temperature, their phase behavior is predictable with only one temperature-independent 
binary interaction parameter (Abu-Eishah, 1991; Moshfegian et al, 1992; Gow% 1993; 
Kleiber, 1994; Lee et ah, 1996; Chimg and Kim, 1997; Lee et ah, 1998). Moreover, it is 
feasible to correlate the binary parameter with temperature for each binary system (Gow, 
1993). Nonetheless, a global correlation applicable for all binaries is difficult to achieve, 
as noted by Abu-Eishah (1991). An alternative approach is to correlate the binary 
parameter with pure component parameters, such as the acentric factor or compressibility 
factor, as shown by Lee and Sun (1992). However, with a growing number of new 
refrigerant mixtures, including alcohol-containing mixtures (Laugier et ah, 1994), and the 
increased ranges of temperature and pressure of interest, this approach is no longer 
adequate for vapor-liquid equilibrium (VLE) information. 

In recent years, much effort has been given to the incorporation of the excess Gibbs free 
energy G*^ models into mixing rule expressions for the attractive term parameter of cubic 
equation of state, as it allows description of vapor-liquid equilibria (VLE) in complex 
systems through a simple EOS (Orbey and Sandler, 1995a;Heidemann, 1996). There is a 
rapidly growing literature exploring techniques for implanting specific G behavior into 
the cubic equation of state through the mixing rules (EOS-G^ framework). The need for 
alternative refrigerant systems brought up the requirement for predictive models over a 
range of conditions, which need only a few experimental data in order to estimate their 
parameters. . Appropriate candidate models that can meet this need are based on the G - 
EOS formalism ( Orbey and Sandler, 1995a). Advantages of this approach are the built-in 
temperatui^ dependence of the parameters of the G^ model, and the ability to utilize 
diverse experimental information, as for example excess enthalpy (H^) experimental data 
(Peng etal, 1995). 

Models that combine equations of state and activity coefficient (or G^) models can be 
categorized into two groups: those that make this link at infinite pressure (Huron and 
Vidal, 1979; Wong and Sandler, 1992), and those that make this link at low or zero 
pressure (Mollerup, 1986; Heidemann and Kokal, 1990; Michelsen, 1990a,b; and Dahl 
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and Michelsen, 1990). Of these, the Wong-Sandler (WS) mixing rules have show 
remarkable extrapolating powers and gained wide acceptance. Wong and Sandler (1992) 
and Wong et al. (1992) incorporated the second virial coefficient as well as solution 
model into the mixing rule obtained by incorporating the excess Helmholtz free energy 
(A'^) of a solution model into that of an EOS at infinite pressure. Available low-pressure 
solution model parameters are found to be useful in this mixing rule because of the 
insensitivity of A^ to pressure. This theoretically correct mixing rule provides a unified 
and consistent way of combining two-parameter cubic equations of state and excess free 
energy models to describe the phase behavior of a wide variety of mixtures and has a 
broad range of applicabiliw. Summarizing in the words of Sandler and Orbey (199S), 
’’Wong and Sandler have developed a mixing rule that combines an EOS with a free 
energy model but produces the desired EOS behavior at both low and high densities 
without being density dependent, uses the existing tables of G“ parameters, allows 
extrapolation over wide ranges of temperature and pressure, and provides a concepmally 
simple method of accurately extending UNIFAC or other low-pressure VLE prediction 
methods to high temperature and pressure”. The original paper has been followed with a 
number of publications by Sandler and his coworkers, including W'ong et al. (1992), 
Pividal et al. (1992), Orbey et al. (1993), Huang et al. (1994), Orbey and Sandler (1995a. 
d), Orbey and Sandler (1996). Orbey and Sandler (1997) and Shiflett and Sandler (199S). 


Orbey and Sandler (1995d) have tested the temperature extrapolation capabilities of the 
conventional vdW and three recently introduced multiparameter mixing rules for a cubic 
equation of state to predict the VLE behavior of some CFG, HCFC, and fluorinated 
alcohol binary mixtures. Tne multiparam^r mixing rules considered were the three- 
parameter mixing rule of Wong and Sandler (1992), a two- parameter empirical mixing 
rule introduced by Stryjek and Vera (1986b) and various other investigators (Adachi and 
Sugie, 1986; Panagiotopoulos and Reid, 1986) and then developed into a multiparameter 
form by Schwartzentruber and Renon (1989), and a new form of a two-parameter 
modified Huron-Vidal mixing rule (Sandler and Orbey, 1995a). Of particular interest was 
determining how well these mixing models could correlate experimental data at a single 
temperature, and then how accurately each could predict phase behavior at other 



temperatures with temperature independent parameters. The results indicate that for some 
new refrigerant mixtures only the excess free energy based equation of state mixing rules 
(Wong and Sandler, 1992, Orbey and Sandler, 1995a) can pro\dde acceptable 
correlations, and also reasonable predictions over a wide range of temperatures. Wliile all 
the models considered performed well for conventional refrigerant mixtures, only the 
excess free-energy based models consistently led to accurate predictions over wide- 
ranges of temperature; the other models result in less accurate predictions and, in some 
cases, false liquid-liquid splitting. 

Recently, loaimidis and Knox (1997) have presented a single energy parameter model 
for several traditionally used and recently proposed refrigerant mixtures based on the 
Wong-Sandler mixing mle. The solution model used is a modification of model based 
on the work of Knox et al. (1984), that assigns one energy and one size related parameter 
per binary mixture. The system studied include fluorocarbons (FCs), CFG, and KCFC 
mixtures and CFC/HCFC-hydrocarbon mixtures. 

Shiflett and Sandler (1998) have modeled seven binary fluorocarbon mixtures forming 
azeotropes using the Wong-Sandler mixing rules with the nonrandom-two liquid (XRTL) 
activity coefficient model (Renon and Prausnitz, 1968). The azeotropes consisted of 
binary mixtures of HFCs, FCs, carbon dioxide, ammonia, and propane. These azeotropes 
demonstrate a variety of nonideal behavior including polar/polar, polar/nonpolar, and 
nonpolar/nonpolar interactions, and in some cases large molecular size differences. Tnese 
systems provided a challenging test for the Wong-Sandler mixing rule. 



CHAPTER 3 


THERMODYNAMIC MODELING 

3.1 Criterion for vapor-Liquid Equilibrium 

The criterion that must be satisfied for components at equilibrium between liquid and 
vapor phases at the same temperature (T) and pressure (P) is given by the equality.' of 
component fugacities throughout the phases. 

fi''=fi'- (3.1) 

The liquid-vapor-phase component fugacities can both be described through fugacity 
coefficients: 

fi^- =Xi 9 i"P = fi'^ = y, 9 >''P (3.2) 


The liquid-phase fugacity coefficient of a component i in a mixture, 9 i'', is given by 

dV-lnZ^ (3.3) 








RT 

v'- 


where Z'^ is the liquid compressibility factor. 

Similarly the vapor-phase fugacity coefficient of a component i in a mixture, 91 ' , is gi\'en 
by 


In^f =-L] 
RT I 


''dP^ 




RT 


W-luZ 


(3.4) 


where is the vapor compressibility factor. 

A suitable pressure-explicit EOS, P=f(V,T) applicable to all pure components and their 
mixtures over the entire density range(vapor to liquid) is required for the solution of Eqs. 
(3.3) and (3.4). 



3.2 Cubic Equation of State 

The PRSV equation of state used in this study is the Peng-Robinson EOS modified by 
Stryjek and Vera(1986a): 
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P = 


RT 


a(T) 


(3.5) 


v-b v{y+b)+b{y-b) 
where v is the molar volume ,R the universal gas constant, a and b the characteristic 
mixture energy and size (covolume) parameters, respectively. Pure component 
parameters a and b are obtained from Eqs. (3.6) through (3.9) from critical temperature 
Tc, critical pressure Pc, acentric factor co, and an adjustable parameter k- characteristic of 
each pure component: 


a{T) = 0.457235 [l + k{l - f 


(3.6) 


where Tr = T/Tc 


k = k,^k,(l + T^%7-T,) (3.7) 

Ato =0.378893 + 1.4897153<2i-0.17131848iy' +0.0196554ni' (3.S) 


RT 

b = 0.077796—^ 
Pc 


(3.9) 


The PRSV equation of state is applicable to a wide 
associating substances. 



of non-polar, polar and 


For the PRSV equation of state and arbitrary set of mixing rules for mixture parameters a 
and b, the fugacity coefficient of component i in a mixture is given by 
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In (p^ = - In 


P{v - b) 
RT 



dnb 
\ J 


Pv 

RT 


+ 


1 

r a 'I 

1 

Iona 

272 ' 


a 

J 


v^a n. j 


v + 6(l — V 2 ) 
+ 6(1 + 72) 


( 3 . 10 ) 


3.3 Mixing Rules: The Wong-Sandler (WS) Model 


The development of Wong-Sandler mixing rules is based on two observations. The first is 
the quadratic composition dependence of the second virial coefficient. This provides one 
of the restrictions on the equation of state a and b parameters, given by 




RTj 


( 3 . 11 ) 


{/ 


For the combining rule, various possibilities exist. In the original %vork, Wong and 
Sandler used 



( 3 . 12 ) 


which introduces a second virial coefficient binary interaction parameter kij (kij=kji). 

The second observation is that the excess Helmholtz free energy of mixing (A^) is much 
less pressure dependent than the excess Gibbs free energy(G ). 

G^[T, .r, , P = low) = {T,x., P = low) =A^{T, x,. , P = high) 

The first of these equalities follows from the fact that G^=A^+PV^, and the PV^ term is 
very small at low pressures. The second of these equalities is a result of the essential 
pressure (or density) independence of A^. The second equation for the a and b parameters 


\ 

( 3 . 13 ) 
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then comes from the condition that 


Alas{T,x.,P = ^)^A^{T,x,,P=^co) = A^{T,x,,P = lo^^)=^G^{T,x,,P = low) (3.14) 

Here the subscript EOS refers to the excess Helmholtz free energy derived from an 
equation of state, while and without subscript refer to a free energy or activity 
coefficient model. Combining these, the basic relation for the Wong-Sandler rule for the 


mixture parameters a and b of a cubic equation of state are 


b-—=yYxi^j b-— 

RT V ; \ RT)y 


(3.11) 


^7 b ^ ‘ b: 


(3.15) 


where is an excess Helmholtz free energy term, and as a result of Eq. (3.14) above, 
any G (activity coefficient) model can be used for A . Here a is a constant dependent on 
the specific EOS selected; (l/2°’^)ln(2°'^-l)=-. 62323 for the PRSV equation used here. 
The parameters in this mixing rule are those of the A^ model and the binary interaction 
parameter ky appearing in the combining rule, Eq. (3.12). 

Equations (3.11) and (3.15), together with Eq. (3.12), completely define a and b in terms 
of A'^oo (the high-density term) and ky (the low-density term). These equations can be 
solved to obtain 




RTcr Y \biRT J 


(3.16) 


RT RTIY ‘ b. cr 


1-D 


(3.17) 
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with Q and D defined as 




^ a ^ 


» J 


RTjy 


\ Ai y 


(3.18) 


and 


Z-J I u D 


a.- At 


■ + ■ 


, b,RT RTo- 


(3.19) 


Using the Wong-Sandler mixing rules given by Eqs. (3.16) and (3.17), the partial 
derivatives of a and b required in Eq. (3.10) for calculating Incp, are: 


and 


d fib 1 

( 1 dn-Q^ 

2 

( d nD'] 

dn, "(l-D) 

[nd ft; J 

{l-Df 

[ 5 J 


1 \ \ d n'^ _ nb ^ ^ 
RT\^ndn.j d n- 

with the partial derivatives of Q and D given by : 


S nD 
5 n, 


1 dn-Q 
n d «; 




RTjy 


d nP _ a-, ^ In/,,. 
5 tti b^RT a 


with 


In V 


1 dnAl 
RTd n. 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 
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3.4 Activity Coefficient Models: Binary Systems 

3.4.1 NRTL (non-random t>vo liquid) Model: 

In Non-Random-Two-Liquid (NRTL) model of Renon and Prausnitz (1968), the molar 
excess Gibbs free energy for a binary mixture is given as 


RT 


= 


Gi^2I ^12^12 




(3.25) 


where 


G2 = 


cr — cr 
i?l2 ^22 


RT 


■ 21 


— oil 

RT 


G,2 =exp(-<2'r,2). 
The activity coefficients are given by 


=exp(-^ 2 rr,,), 


In/i =xl 



f Gix 

'I 

2 

1 ^ 12^12 

^21 

^ 21 * 



^ 2 V 

('^2 ^12 “^1 } 


(3.26) 


In^i =^i" 


' G.2 




V'^2 +G:,2-Ly 




Gi^21 


{x^ +G2,X,)' 


(3.27) 


It is a three parameter (a, xu, and xii) model in which a is a non-randomness parameter. 


3.4.2 UNIQUAC (UNIversal QUAsi Chemical) Model: 

The UNIQUAC model of Abrams and Prausnitz (1975) is 

qE qE G^ 

= (combinatorial) ^ (residual) 

RT RT^ ^RT 


(3.28) 


where 
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RT 


{combinatorial) = Xj In — + In — + — 

.r, X 2 2 


1 1 

^i-tj In — + ^,'■^2 

V n Yi 


(3.29) 




i?r 


{residual) = ln(^, + 9^ r,] ) - ^ 2'''^2 '^v ) 


(j)i = segment or volume fraction of component i 


f \ 

XQ ■ 

0i = area fraction of component i = ■ - — 

L^j^j 
V J J 

P =volume parameter of component i; 
qj = surface area parameter of component i; 


f \ 

V i J 


( a.A 


V ^ J 


= adjustable parameter; 


(3.30) 


r.=exp{-!;^} = exp| 

Uij = average interaction energy for the interaction of molecules of component i with the 
molecules of component); and 
z = coordination number which is usually taken as 10. 

The activity coefficients are given by 


In/; = In/f {combinatorial) -\-\n /f {residual) 


(3.31) 


A'. 2 (p. ^ J J 


In/f = q, 


( \ 


-X, 






k ^kj 


(3.32) 


(3.33) 
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where /. =L{r. - q;)-{r.-l) 

The structural parameters r; and qi can be calculated as the sum of the group volume and 
area parameters Rk and Qk. 

(3.34) 

k 

(3.35) 


where is the number of groups of type k in a molecule of component i. 



CHAPTER 4 


RESULTS AND DISCUSSION 


4.1 Pure Component Constants and Mixture VLE Data 

Table 4.1 gives the literature source of isothermal VXE data for 37 binary refrigerant 
mixtures. The critical temperature (Tc), critical pressure (Pc), and acentric factor (co), for 
each component in the mixture, taken from the references given in Table 4.1, are 
provided in Table 4.2. A computer program was used to evaluate the ki constant of the 
PRSV equation of state by minimizing the overall deNiation in vapor pressure (AAD). 
The ki constant and the AAD are both provided in Table 4.2 together with the literature 
source of vapor pressure data. The overall deviation in pressure is calculated using the 
following formula 


aad=-Y 


where N is the number of data points, Pexp 
pressure. 


ioo (4,1) 

p 

exp 

is measured pressure, and Pcai is the calculated 


4.2 Data Reduction 


In the Wong-Sandler mixing rule the EOS parameters for a homogeneous liquid or vapor 
mixture is computed from Eqs. (3.16), (3.17) together with Eq. (3.12). The single 
adjustable parameter ky [Eq. (3.12)] for each binary pair is referred to as the WS mixing 
rule parameter (ki 2 ). Three excess free energy (or equivalently activity coefficient) 
models were evaluated for calculating G^: NRTL, UNTQUAC, UNIQUAC-R. Parameters 
ti 2 ,T 2 i for each binary mixture are referred to as the model parameters. The values of 
the pure component parameters q, qi in the UNIQUAC model were taken from loannidis 
and Knox (1997) or calculated from the group contribution parameters R^, Qk, 
respectively, given by Fredenslund and Sorensen (1994). The NRTL model has the 
additional advantage of having a third adjustable parameter, a. This is typically set to a 
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Table 4.1 Literature source of experimental VLE data for refrigerant mixtures (contd.). 


No 

System 

N 

T(K) 

P Range (MPa) 

Reference 

1. 

R-32+R-143a 

13 

256.00 

0.35-0.45 

Shiflert and Sandl 






(1998) 

2. 

R-32+R-116 

14 

253.55 

1.056-0.4108 

Shiflert and Sandl 






(1998) 

3. 

R-32+R-125 

11 

257.85 

0.4008-0.4815 

Shiflert and Sandl 



14 

284.40 

0.9459-1.1477 

(1998) 

4. 

R-32+R-290 

11 

257.85 

0.2863-0.4815 

Shiflert and Sandl 






(1998) 

5. 

R-23+R-116 

12 

209.60 

0.2094-0.2681 

Shiflert and Sandl 



15 

227.60 

0.4398-0.5737 

(1998) 

6 . 

R-32+R-134a 

5 

263.15 

0.1997-0.579 

Chung and Kim 



5 

273.15 

0.2916-0.810 

(1997) 



5 

283.15 

0.414-1.104 




5 

293.15 

0.5709-1.471 




5 

303.15 

0.769-1.9227 




5 

313.15 

1.015-2.474 




5 

323.15 

1.316-3.136 


7. 

R-125+R-134a 

6 

268.15 

0.254-0.544 

Kobayashi and 



8 

273.15 

0.308-0.633 

Nishiumi (1998) 



4 

278.15 

0.415-0.689 




8 

283.15 

0.442-0.882 




7 

293.15 

0.606-10158 


8 

R32+R-125 

6 

283.05 

0.915-1.104 

Chung and Kim 



5 

293.05 

1.202-1.48 

(1997) 



5 

303.05 

1.563-1.923 




4 

313.05 

2.001-2.477 


9. 

R-22+R-12 

11 

356.15 

2.413-3.88 

Nishumi et al.(1995', 



18 

343.81 

1.869-3.01 


10. 

R-22+R-123 

13 

414.57 

1.785-5.813 

Nishumi etal.(1995] 



13 

383.15 

0.985-5.158 




14 

353.92 

0.49-3.216 




12 

333.70 

0.279-2.381 


11. 

R-134a+R-124 

6 

296.45 

3.7-6.21 

Lee et al. (1996) 



9 

302.25 

4.31-7.48 




8 

307.25 

4.98-8.65 

1 

12. 

R-124+R-142b 

12 

298.15 

3.38-3.79 

Lee et al. (1996) 



11 

312.15 

5.08-5.73 



Table 4.1 Literature source of experimental VLE Data for refrigerant mixtures (contd.). 


No 

System 

N data 

T(K) 

P Range (MPa) 

Reference 

13. 

R-142b+R-141b 

13 

323.4 

0.185-0.6846 

Kang and Lee (199t 

14. 

R-142b+R-140a 

13 

323.4 

0.0465-0.6846 

Kang and Lee (199( 

15. 

R-141b+R-140a 

14 

323.4 

0.0465-0.185 

Kang and Lee (199( 

16. 

R-32+R-22 

13 

283 

0.6722-1.1035 

Kang and Chung. 






(1996) 

17. 

R-32+R-12 

16 

283 

.4204-1.1035 

Kang and Chung 






(1996) 

18. 

R-32+R-40 

13 

283 

.3619-1.1035 

Kang and Chung 






(1996) 

19. 

R-134a+R-141b 

5 

278.15 

0.101-0.269 

Zheng et al.(1990) 



6 

288.15 

0.101-0.389 




7 

318.15 

0.221-0.999 




11 

333.15 

0.34-1.504 


20. 

R-23+R-13 

16 

272.944 ■ 

1.96-2.49 

Stein and Proust. 



18 

249.66 

1.218-1.499 

(19-1) 



18 

224.611 

.446-.513 




21 

199.6388 

.154-. 164 


21. 

R-134a+R-123 

6 

289.17 

0.12-0.393 

Kubota et al.(1991) 



7 

303.15 

0.201-0.628 




8 

318.15 

0.278-1.096 




10 

332.74 

0.476-1.495 


22. 

R-22+R-114 

5 

253.15 

0.076-0.229 

Kubota et al.(1990a: 


i 

8 

283.15 

0.181-0.619 




12 

313.15 

0.616-1.175 




19 

338.15 

0.799-2.594 


23. 

R-12+R-114 

10 

283.15 

0.18-0.394 

Kubota et al.(1990a) 



10 

313.15 

0.387-0.896 


24. 

R-13B1+R-12 

5 

253.15 

0.207-0.411 

Kubota et al. (1990a] 



18 

283.15 

0.463-1.066 




12 

313.15 

1.069-2.087 


25. 

R-22+R-142b 

15 

263.15 

0.098-0.355 

Kubota et al. (1990b) 



16 

273.15 

0.146-0.498 




16 

283.15 

0.21-0.679 




21 

293.15 

0.293-0.824 




18 

313.15 

0.526-1.534 




13 

338.15 

0.995-2.704 


26. 

R-32+R-142b 

8 

295.45 

0.31-1.571 

Lee et al.(1998) 



8 

304.55 

0.409-1.998 




8 

314.95 

0.549-2.558 
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Table 4.1 Literature source of experimental VLE data for refrigerant Mixtures. 


No 

System 

N data 

T(K) 

P Range (MPa) 

Reference 

27. 

R-32+R123 

9 

294.95 

.0816-1.554 

Lee et al.(1998) 



9 

304.55 

.116-1.998 




10 

313.95 

.159-2.527 


28 

R-134a+R600a 1 

17 

293.66 

0.3067-0.58 

Bobbo et al. (1998 


1 

11 

303.68 

0.41-0.7809 


29. 

R-134a+R-236fa i 

9 

283.62 

0.162-0.42 

Bobbo et al. (1998 



10 

303.68 

0.325-0.78 


30. 

R600a+R236fa 

15 1 

303.68 

0.325-0.41 

Bobbo et al. (1998' 

31. 

R-22+R23 

8 

273.15 

0.7152.355 

Roth et al. (1992) 



7 

293.15 

1.11-3.85 




8 

303.15 

1.525-4.8 




8 

323.15 

2.4-5. 2 


32. 

R-152a+R-150a 

10 

303.2 

.134-.685 

Lim et al. (1996) 



10 

323.2 

.158-1.147 


33. 

R-134a+R-142b 

i 10 

268.0 

.1189-.2425 

Klieber. (1994) 



i 11 

283.0 

.2066-.4132 




' 11 

298.0 

.3365-.6622 


34. 

R-134a+R-12 

1 10 

258 

.1815-.1631 

Klieber (1994) 



10 

278 

.3602-.3485 




10 

298 

.6477-.6622 


35. 

R-134a+R-152a 

1 8 

255 

.1306-.1438 

Baieber(1994) 



I 8 

275 

.2821-.3129 




1 9 

298 

.5936-.6622 


36. 

R-134a+R-116 

1 10 

251 

.9804-.1210 

Klieber (1994) 



i 8 

275 

1.1709-.3129 


37. 

R-14+R-13 

1 12 

199.80 

.1544-1.5272 

Proust and Stein. 



1 



(1979) 

38. 

R-22+R-134a 

i 7 

273.16 

.2925-.4983 

Arita et al. (1991) 



i 6 

298.16 

.6651-1.0445 




1 7 

323.16 

1.3173-1.9420 





4.2: Critical temperatures and pressures, acentric factor, ki and overall pressure deviation (contd.). 
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fixed value of 0.3 for vapor-liquid equilibria. In the case of few binary mixtures such as 
R-32/R-13'4a, R-32/R-125 and R-32/R-290, an a value of 0.2 was found to pro\ide a 
better fit which is consistent with the highly nonideal liquid-liquid equilibria at high 
temperatures. A three parameter fit of the isothermal \nLE data using the WS mixing 
rules with the models (k^, tn, X 21 ) was performed. The computer program used a 
simplex algorithm to optimize model parameters (ti; and X 21 ) and the WS mi.xing rule 
parameter (ki 2 ). The VLE calculations presented here were done using an isothermal 
bubble point algorithm. A flow diagram of the algorithm is presented in Fig. 4.1. When 
optimizing the model parameters, the following objective function was used 

F = y|p. -p. .! (4.2) 

i 

for the minimization in a simplex formalism. 

The WS model showed good convergence ability, which is of vital importance in VLE 
calculations, due to its robusmess to parameter variation. The program used for 
correlating the data required an initial guess for each parameter. The initial values of 
ti 2 , T 21 and kia were chosen to be 0.1. The program was run and fitted parameters along . 
with the calculated pressure and vapor composition were computed. To ensure that a 
global vs. a local minimum convergence was obtained, values of X 12 , X 21 and k;: were 
varied in increments of 0.01. The calculated X 12 , xii and ki 2 which provided the best fit 
with the experimental data were identified as the optimum set of parameters and are 
presented in Table 4.3. The optimum set of UNTQUAC parameters were used also in 
VLE calculations with the UNIQUAC-R model to assess the sensitivity of the entropic 
term (combinatonal part) on phase equilibrium predictions. 

4.3 Results and Discussion 

Tables 4.4 through 4.8 give the VLE results for five systems for which only P-x data 
were available; measured and calculated pressures and the calculated vapor compositions 
are presented along with percent average absolute deviation in pressure. Figures 4.2 
























Table 4.3: Calculated EOS, NRTL and UNIQUAC parameters (contd.). 
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Table 4.5: Vapor-liquid equilibrium pressures and phase compositions for R-32{1)/R-116(2) mixtures at 
with different models. 
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able 4.6:Vapor-liquid equilibrium procures and phase compositions for R-32{1)/R-125{2) mixtures 


o 

“O 

o 


E 


LU 

O) 

c 

0) 

I 

SZ 




ir 

1 

o 

< 

3 

a 

z: 

3 


0.942 

0.982 

1.016 

1.044 

1.067 

1.086 

O 1 

o • 

1.112 

1.128 

1.138 

1.143 

CO 1 

t — ' 

T” 

r— 

1.148 

o 

CD 

d 

P(MPa) 

o 

03 

o 

O 

< 

3 

O 

z: 

3 


0.942 

CM 

CO 

CD 

d 

1.016 

1.044 

1.067 

1.085 

o 

o 

1.111 

1.128 

1.137 

1.143 

1.146 

1.147 

1.148 

0.181 



NRTL 


0.942 

0.984 

1.017 

O 

— 

1.067 

1.084 

1.099 

1.111 

1.128 

1.139 

1.145 

1.148 

1.150 

1.148 

CO 

CO 

o 

d 


a. 

X 

CD 



0.946 

9860 

1.018 

1.044 

1.066 

CO 

CD 

O 

1.098 

1.110 

CO 

CM 

1.139 

CD 

1.149 

o 

ur> 

1.148 

Q 



UNIQUAC-R 

CO 

CM 

I 

o 

o 

o 

di 

0.137 

0.250 

CO 

d 

0.423 

0.489 

0.545 

9690 

0.681 

0.753 

0.818 

0.877 

0860 

o 

o 

o 

x— 

i 


calc 

UNIQUAC 


o 

o 

o 

d 

0.136 

0.250 

CO 

CO 

d 

0,422 

CO 

OD 

d 

0.545 

9690 

0,681 

0.753 

0.818 

0.877 

0.930 

1.000 




NRTL 


o 

o 

o 

d 

0.136 

0.247 

0.339 

0.416 

0.483 

0.540 

0.592 

0.681 

0.754 

0.819 

0.877 

CD 

CN 

CD 

C 

1.000 


x" 


c 

c 

c 

c 

1 CO 

> c 

) T- 

ijC 

0 204 

0 289 

^ cc 
1 cc 

( CO 

! c 

0 435 

1 

> CD 
■ ^ 

i C 

. ^ 

> u: 

■ IT, 

> c 

' 

) IT. 
) <X 

> c 

' OC 
) CO 
) N 

) C 

) x- 

. s 

> c 

. ^ 

) CC 

) c 

' a 

. D 

> a 

i C 

) C 
4 C 
) C 

) X- 

) 

I 

) 


oAAD = average absolute deviation- defined a.v(^a'/>.v((exp-ca/c)/e>q3)*10o)/ 




Table 4.6: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-125{2) mixtures 
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fable 4.8: Vapor-liquid equilibrium pi^ssures and phase compositions for R-23{1)/R-116{2) mixtures at 209.6 K 
with different models (contd.). 
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able 4.8: Vapor-liquid equilibrium pressures and phase compositions for R-23(1)/R-116(2) mixtures 
with different models. 
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Table 4.9: Vapor-liquid equilibrium pressures and phase compositions for R-32{1)/R-134a(2) mixtures 
with different models (contd.). 
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Table 4.9: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-134a(2) mixtures 
with different models (contd.). 
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Table 4.9: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-134a(2) mixtures 
with different models. 
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Table 4.10: Vapor-liquid equilibrium pressures and phase compositions for R-125(1)/R-134a(2) mixtures 
with different models (contd.). 
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fable 4.10: Vapor-liquid equilibrium pressures and phase compositions for R-125{1)/R-134a(2) mixtures 
with different models. 


P(MPa) 


UNIQUAC-R 1 


0.581 

9090 

6990 

0./32 

0.810 

0.901 

0.983 

1.158 

1.198 

0.150 1 

0|B0 

UNIQUAC 


0.581 

CD 

O 

CD 

d 

6990 

0.732 

0.810 

0.901 

0.983 

1.158 

1.196 

0.150 


NRTL 

293.15 K 

0.581 

0.604 

9990 

CM 

CO 

d 

0.813 

9060 

6860 

o 

CD 

T — 

1.196 

CO 

O 

CO 

d 


expti 


0.581 

CD 

O 

CD 

d 

8990 

0.732 

0.816 

8060 

€860 

1.158 

1.196 

%AAD 

\ 

calc 

UNIQUAC-R 

o 

o 

o 

d 

O 

d 

CD 

CO 

CN 

d 

■ 0.394 

0.535 

0.671 

0.775 

8960 

o 

o 

o 

1.9446 

UNIQUAC 

o 

o 

o 

d 

o 

d 

0.236 

0.394 

0.535 

0.671 

0.775 

896 0 

o 

o 

o 

— 

1.945 

NRTL 

o 

o 

o 

d 

0.101 

0.235 

0.398 

0.540 

CD 

d 

0.776 

0,957 

o 

o 

o 

1.916 

expti 


o 

o 

o 

d 

0 097 

0.230 

0.386 

0 535 

0.662 

OZZO 

8960 

000 T 

%AAD 

X 


o 

c 

c 

c 

LD 

ur. 

c 

c 

0.137 

0 255 

cc 

<r 

c 

0.526 

h- 

IC 

CD 

C 

0 931 

1000 



VoAA D = average ahsohite deviation — defined as <7/).v((exp calc)! exp ) * 1 Oo)/ N 







^ i 




CN 

1C 

CM 




CM 

CO 

I 

tr 


o 

CO 

c 

o 


CO 


o 


Q. 


£ 


o 


o 


CD 


CO 


CO 


SI 


Cl 


ID 


c 


CO 


CO 


o 


k« 


D 


(0 

•d 

CO 

'4-* 

p 

c 

Lm 

Q. 

o 

E 

3 

CO 

*JL- 

o 


"O 


o 

3 

£ 

cr lu 

o 

D) 


C 

'5 

o 


w 

o 

"T 

v« 


O 

*"D 

a 

CO 

JZ 


I 


o 

n 

«3 


P(MPa) 


UNIQUAC-R 


0.904 

1.018 

1.058 

1.080 

1.090 

1.098 

1 .409 1 

O 

CO 

o 

UNIQUAC 


0.904 

1.017 

1.058 

6Z0T 

1.090 

CD 

CD 

O 

T — 

1.410 

NRTL 


0.904 

8660 

1.058 

1.098 

m 

o 

T— 

1.098 

CD 

00 

to 

d 


expti 


lO 

o 

cd 

1 CO 
CM 

! 

0.905 

0.992 

1.058 

1.102 

1.125 

1.104 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

o 

CD 

00 

CM 

d 

00 

d 

ozro 

0.900 

o 

o 

o 

1.106 

NRTL UNIQUAC 

o 

o 

o 

CD 

0.285 

CO 

CO 

d 

0.719 

CD 

CD 

CO 

d 

o 

o 

q 

\ 

^ 1.107 

o 

o 

o 

d 

0.292 

0.503 

0.725 

CO 

CD 

CO 

d 

1.000 

0.923 

expti 


o 

o 

o 

o 

0.283 

009 0 

0.735 

Z680 

1.000 

aw% 



c: 

c 

c 

c 

0.241 

0.457 

0.710 

0 895 

O 

O 

C 



1.201 

1.345 

1.396 

1.427 

1.448 

1.464 

1.232 1 

1.201 

q 

x— 

in 

CD 

CO 

CM 

CO 

tT 

xr 

CD 

'M' 

CO 

CM 

1.201 

CD 

T— 

q 

LO !CO 
CDilO 

T— ; ■rr* 

00 

CD 

(D 

T 

0.378 

1.202 

o 

CO 

T“ 

UD f CD 
CDiLO 

q^qr 

i 

o 

00 

CM 

c-. 

%AAD 

0.000 

0.280 

1 

} 

CCMt- 
N-.CM 
'O” 11^ 

c id 

CO 

o 

CD 

d 

o 

o 

q 

0.922 

O 

o 

o 

d 

0.279 

r^ 

r^ 

xr 

o 

i 

5 cm 
jd 

8060 

o 

o 

q 

T“ 

CM 

CD 

d 

o 

o 

o 

d 

0.286 

r^ 

CD 

C 

jcD 

jCM 

id 

i 

i 

968 0 

o 

o 

q 

0.702 

o 

o 

o 

d 

0 279 

i 

COiCM 
CD ICO 

did 

' 1 

Z680 

o 

o 

q 

aw% 

o 

o 

c 

c 

T" 

o 

c 

i 

' CO JCD 
" UO ■ CD 

) die 

n 896 

‘ C 

1 C 

1 c 



Vo A AD — average ahsoluie deviation — defined as ahs{(Qxp c*a/c)/exp)* lOOj/ A/ 


bfe 4.11: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-125(2) mixtures 
with different models. 
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Fig. 4.10 Vapor-liquid equilibrium pressures and phase compositions for 
R-32(1 )/R-1 25{2) mixtures at different temperatures (K). 
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Table 4.12: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-12(2) mixtures 
with different models. 




UNIQUAC-R 


1.865 

CO 

lO 

CT> 

T— 

2.050 

CO 

04 

2.257 

5 

CO 

cvi 

2.392 

2.471 

2.588 

2.668 

2.797 

2.826 

2.855 

2.874 

2.937 

2.979 

2.999 

3.045 

0.472 1 

P(MPa) 

calc 

UNIQUAC 


1.877 

1.953 

2.050 

XT 

CO 

X— 

CN 

2.257 

2.341 

2.392 

2.471 

2.588 

00 

CD 

CD 

CN 

2,797 

2.826 

2.855 

2.874 

2.937 

2.979 

2.999 

3.045 

0,483 



NRTL 


1.877 

1.953 

2.050 

2.184 

2.257 

2.341 

2.392 

2.471 

2.588 

2,668 

2.797 

2.826 

2.855 

2.874 

2.937 

2,979 

2.999 

3.045 

CO 

00 

d 


expti 



cn 

CD 

00 

1.947 

2.026 

CO 

CN 

2.252 

2.341 

2.405 

2.479 

2.612 

2.675 

2.798 

2.818 

2.853 

2.872 

2.916 

2.951 

2.970 

o 

o 

CO 

%AAD 



UNIQUAC-R 

343.81 K 

o 

o 

o 

d 

0.059 

00 

CO 

d 

0.236 

CD 

CN 

d 

0.359 

o 

o 

d 

CD 

d 

0.562 

0.631 

0.748 

0.775 

CO 

o 

CO 

d 

0.821 

00 

oo 

d 

0.928 

0960 

o 

o 

o 

1.827 


calc 

UNIQUAC 


o 

o 

o 

d 

6900 

CO 

CO 

d 

CD 

CO 

CN 

d 

0.293 

0.359 

o 

o 

d 

0.464 

0.562 

0.631 

0.748 

9ZZ0 

CO 

o 

00 

d 

0.821 

00 

CO 

d 

00 

CN 

a> 

d 

0960 

1.000 

1.726 



NRTL 


o 

o 

o 

d 

6900 

CO 

CO 

d 

0.236 

0.293 

0.359 

o 

o 

d 

CD 

d 

0.562 

CO 

CD 

d 

0.748 

0.775 

0.803 

0.821 

00 

00 

d 

00 

CN 

CD 

d 

0960 

o 

o 

o 

CD 

00 

CD 

CN 


CL 

X 

0 


o 

o 

o 

d 

0.057 

0.132 

0.243 

0.301 

0.371 

0.414 

0.481 

0.574 

CO 

00 

CO 

d 

0.750 

0.776 

0.802 

0.820 

0.881 

0.925 

0.947 

1.000 

%AAD 

} 


o 

c 

c 

c 

0 036 

cc 

oc 

c 

c 

c\ 

T“ 

c 

CC 

c\ 

c\ 

c 

CC 

CC 

c\ 

c 

0.327 

0 393 

0498 

0 575 

0.708 

0.739 

c\ 

N- 

h- 

C 

CO 

a: 

c 

999 U 

0 917 

0.942 ^ 

c 

c 

c 



%AAD — average absolute deviation — defined as abs{{exp calc ) / exp ) 


m 

0 

3 

' 4 -* 

X 

1 

CM 

T~ 

(k 


CM 

CM 

I 

a: 


m 

c 

o 


(/> 


o 


a 


E 


o 


o 


o 


</) 


ns 


JC 


a 


TS 


c 


ns 


CO 


o 




3 

>0^ 

CO 

CO 

xs 

CD 

c 

u. 

Q. 

o 

-S- 

S 

3 

i£ 

*3 

CD 

£ 

xs 


O 

3 

E 

cr 

LU 

CD 

O) 

“O 

•4>«( 

c 

'5 

o 




CD 

"T 

te 

O 

3 

a 

ns 

> 

JZ 

4»<ji 

’i 



'T* 




CD 


21 







P(MPa) 


UNIQUAC-R 


1.784 

1.981 

2.238 

2.814 

3.022 

3.226 

3.387 

3.810 

4.001 

4.458 

4.649 

4.889 

5.089 

r-- 

GO 

CD 

o 

CD 

O 

UNIQUAC 


1.795 

1.991 

2.243 

o 

(N 

00 

csi 

3.024 

3.224 

3.388 

3.830 

CD 

O 

O 

4.453 

CD 

4.893 

5.091 

2.205 


NRTL 


1.795 

1.949 

2.164 

2.713 

2.922 

3.131 

3.305 

3.784 

3.978 

4.472 

4.652 

5.018 

5.207 

8600 


exptl 


414.57 K 

1.785 

CO 

CD 

CO 

2.158 

2.687 

2.922 

3.104 

CO 

CO 

CO 

3.751 

CD 

O 

O 

4.501 

4.653 

4.884 

CO 

00 

T— 

to 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

0.091 

CD 

00 

d 

CO 

CO 

d 

0.391 

0.428 

0.456 

0.520 

0.542 

8890 

0.601 

CD 

d 

0.620 

9.016 

UNIQUAC 

o 

o 

o 

d 

0.091 

0.187 

0.347 

0.391 

0.429 

0.457 

0.520 

0.542 

8880 

o 

CD 

d 

0.645 

0.630 

9.2 

NRTL 


o 

o 

o 

d 

O 

CO 

o 

d 

0.175 

0.350 

o 

o 

d 

CO 

d 

0.475 

0.547 

0.571 

0.622 

0.637 

00 

CD 

CD 

d 

0.673 

10.471 

exptl 


o 

o 

o 

o 

0.053 

CO 

CD 

d 

0.325 

0.375 

0.411 

0.458 

0.501 

0.533 

0.583 

0.593 

0.599 

869 0 

%AAD 

X 


c 

c 

c 

c 

0 034 

cr: 

h- 

c 

c 

cc 

X“ 

c 

T~ 

CN 

Cs 

c 

n 957 

0.286 

CO 

CD 

CO 

c 

CO 

cr 

lo 

c 

0468 

LC 

cr 

c 

c 
ir 
- ID 

) C 

<7 

N 

IT 

C 



%AA D ~ average absolute deviation ~~ defined as a/).v((exp - calc) / exp) * 1 0O)/ hi 



Table 4.13: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-123(2) mixtures 
with different models (contd.). 
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Table 4.13: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-123(2) mixtures 
with different models (contd.). 
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Table 4.14: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-124{2) mixtures 
with different models (contd.). 
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Table 4.15:Vapor-liquid equilibrium pressures and phase compositions for R-124a(1)/R-142b(2) mixtures 
with different models (contd.). 
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Table 4.17: Vapor-liquid equilibrium pressures and phase compositions for R-142b(1)/R-140a{2) mixtures at 323.4 K 
with different models. 
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4.18: Vapor-liquid equilibrium pressures and phase compositions for R-141b(1)/R-140a(2) mixtures at 323.4 K 
with different models. 
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able 4.19: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-22(2) mixtures at 283.0 K 
with different models. 


P(MPa) 


UNIQUAC-R 

0.675 

0.731 

0.777 

0.821 

6990 

0.910 

0.950 

0.988 

0.911 

0.988 

1.024 

1.062 

1.101 

Dsi 

CM 

d 

0|B0 

UNIQUAC 

0.675 

0.731 

ZZZO 

0.821 

698 0 

0.910 

0.950 

CO 

CO 

CD 

d 

0.911 

8860 

1.024 

1.062 

1.101 

oo 

CM 

CM 

d 


NRTL 

0.675 

0.729 

0.774 

0.817 

998 0 

0.908 

0960 

0660 

6060 

0660 

CO 

CM 

O 

1.065 

1.101 

0.206 


expti 


0.677 

0.733 

0.777 

o 

CM 

CO 

d 

<D 

CO 

d 

CO 

o 

CD 

d 

0960 

0660 

0.904 

6860 

1.026 

1.065 

1.104 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

0.173 

0.299 

CO 

o 

d 

0.520 

0.610 

CD 

CD 

d 

0.775 

0.611 

0.775 

0.851 

0.928 

o 

o 

o 

1.089 

UNIQUAC 

o 

o 

o 

d 

0.173 

0.299 

CO 

d 

0.520 

0.610 

9690 

9ZZ0 

0.612 

9ZZ0 

0.851 

0.928 

o 

o 

o 

1.089 

NRTL 

o 

o 

o 

d 

0.171 

0.295 

CO 

o 

d 

0.516" 

CO 

o 

CD 

d 

0.694 

9ZZ0 

CD 

O 

CD 

d 

9ZZ0 

0.851 

0.927 

o 

o 

o 

0.658 n 

expti 

000 0 

0.168 

0.292 

0.399 

lO 

d 

0.604 

0.692 

0.772 

9090 

0.787 

0.851 

0.928 

1.000 

%AAD 

>r 

O 

C 

C 

c 

T“ 

c 

n 208 

CD 

CD 

C\ 

c 

0 406 

0.503 

c\ 

c 

C£ 

C 

0 703 

c 

in 

c 

c\ 

c 

D- 

c 

c 

C 

CC 

c 

0 901 

1 000 



Vo/l/i D = iivc'raj^e ahsohUe devicU ion — defined as {^ahs{(exp calc)/ exp)* ]0Q)/ N 



Table 4.20: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-12(2) mixtures at 283.0 K 
with different models. 
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Table 4.21: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-40{2) mixtures at 283.0 K 
with different models. 


P(MPa) 


UNIQUAC-R 

0.362 

0090 

9090 

0.704 

8ZZ0 

0.843 

0.901 

0.954 



0.844 

0.954 

1.006 

1.054 

1.100 

0.190 1 

calc 

UNIQUAG 

0.362 

009 0 

0.605 

0.704 

0.778 

0.844 

0.901 

0.955 

0.845 

9960 

1.006 

1.054 

1.100 

0.191 


1 — 
tr 
z 

0.362 

0.494 

6690 

OOZO 

0.776 

CO 

d 

0.903 

Z960 

0.845 

896 0 

1.009 

1.056 

O 

O 

T — 

0.405 


expti 


0.362 

909 0 

ID 

O 

CD 

d 

0.704 

ZZZO 

0.844 

0.901 

9960 

0.841 

9960 

00 

o 

o 

1.055 

1.104 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

0.327 

CO 

d 

0.598 

0Z90 

o 

CO 

d 

0.782 

0.831 

CO 

d 

0.832 

0.883 

0.938 

o 

o 

o 

'T- 

co 

CM 

d 

UNIQUAC 

o 

o 

o 

p 

'\ 

N- 

CN 

CO 

d 

\ 

00 

d 

9690 

0.671 

0.730 

0.782 

CM 

CO 

CO 

d 

0.731 

0.832 

CO 

00 

00 

d 

0.938 

o 

o 

o 

0.278 

NRTL 

o 

o 

o 

d 

0.319 

0.479 

9690 

6990 

o 

CO 

d 

0.782 

0.831 

CO 

d 

0.831 

CO 

00 

CO 

d 

0.937 

o 

o 

o 

T — 

0.319 

expti 

0000 

0.324 

0.481 

Z690 

6990 

0.731 

98Z0 

0.834 

0.731 

0.834 

Z880 

0 940 

1.000 

%AAD 

x" 

O 

O 

O 

C 

n 098 

0.191 

c 

c 

CO 

c 

■ 0 399 

009 0 

c 

c 

CD 

C 

0.700 

c\ 

c 

ic 

c 

OOZO 

c\ 

c 

oc 

c 

'T- 

c 

cr 

C 

000 t 



%AAD = Cfvef'oge absolute deviation — defined as (^^cibs(^GKp calc) /Q'xp)* 100^/ N 



Table 4.22: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)R-141b{2) mixtures 
with different models (contd.). 

— — 1 

P(MPa) 


UNIQUAC-R 


0.035 

0.109 

0.131 

0.165 

0.212 

0.270 

0.348 

CD 

CO 


0.053 

o 

r— 

d 

S) 

r” 

d 

rsi 

d 

d 

0.308 

0.38^ 

0.485 

0.452 i 

calc 

UNIQUAC 


0.035 

0.110 

0.132 



0.166 

0.213 

0.271 

0.350 

2.080 


0.053 

0.102 

0.154 

0.211 

0.270 

0.309 

0.388 

0.484 

CO 

CO 

CD 

d 


NRTL 


0.035 

0.101 

0.128 

0.165 

_____ 

0.219 

0.271 

00 

CO 

d 

1,790 


0.053 

CD 

CO 

O 

d 

0.132 

CD 

d 

0.270 

0.318 

CD 

CD 

CO 

d 

0.489 

5.403 


expti 


278.15 K 

0.035 

0.101 

0.132 

0.167 

0.213 

0.269 

00 

CO 

d 

%AAD 

288.15 K 

0.053 

0.101 

0.154 

0.215 

0.270 

CD 

O 

CO 

d 

0.389 

CD 

CO 

d 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

I 0.737 

0.792 

0980 

0.901 

0.945 

1.000 

4.501 

o 

o 

o 

d 

0.518 

0.702 

0.803 

CD 

CO 

d 

0.894 

0.942 

o 

o 

o 

4.515 

UNIQUAC 

o 

o 

o 

d 

CO 

d 

0.792 

0.851 

0.902 

0.945 

O 

O 

O 

5,425 

o 

o 

o 

d 

0.518 

0.702 

CO 

o 

CO 

d 

to 

CD 

00 

d 

9680 

0.943 

1.000 

4.548 

NRTL 

o 

o 

d 

0.722 

0.791 

6980 

0.912 

0.949 

O 

O 

q 

5.469 


o 

o 

d 

0.432 

0.665 

o 

CO 

d 

0880 

0.911 

0960 

O 

O 

q 

5.894 

expti 


o 

o 

o 

d 

0.647 

0.759 

0.819 

ZZ8 0 

0.925 

1.000 

avv% 


o 

o 

q 

d 

0.498 

999 0 

0.768 

0 834 

Z98 0 

v-" 

s 

o 

o 

o 

o 

%AAD 

X 



o 
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c 
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c 
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c 

C\ 

IC 

•r- 
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C\ 

c 

CC 
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i C 

c\ 

c: 
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j c 

cx: 

I T- 

-Jn 

)IC 
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c 



%AAD = average ahsolule deviation - defined a.';i^abs{{ey:p- calc)l e^)*\00)l N 



Table 4.22: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)R-141b(2) mixtures 
with different models (contd.). 

— 1 

P(MPa) 


UNIQUAC-R 


0.164 

0.217 

0.335 

0.431 

0.533 

0.637 

0.798 

0.980 

1.155 

CD 

D 

CD 

d 

0|B0 

UNIQUAC 


0.154 

0.217 

LO 

CO 

CO 

d 

0.432 

0.534 

CO 

CO 

CD 

d 

66Z0 

0.981 

1.156 

o 

CD 

d 

NRTL 


0.154 

0.215 

CO 

CO 

CO 

d 

0.431 

0.535 

CX) 

CO 

CD 

d 

66Z'0 

0.981 

1.157 

0.781 


expti 


318.15 K 

0.156 

0.221 

CD 

CO 

CO 

d 

0.431 

0.534 

00 

CO 

CD 

d 

o 

o 

00 

d 

6660 

1.158 

%AAD 

>; 

calc 

UNIQUAC-R 

o 

o 

o 

d 

o 

o 

CO 

d 

0.571 

CD 

CO 

CD 

d 

89Z0 

CM 

CM 

00 

d 

Z880 

0.945 

o 

o 

o 

11.536 

UNIQUAC 

000 0 

o 

o 

CO 

d 

0.571 

CD 

CO 

CD 

d 

0.766 

CM 

CM 

00 

d 

OO 

CO 

d 

CD 

CX) 

d 

o 

o 

o 

11.595 

NRTL 

O 

O 

O 

d 

0.295 

893 0 

0.684 

CD 

d 

0.820 

0.884 

0.944 

o 

o 

o 

10.948 

expti 


o 

o 

o 

d 

0.215 

0.604 

0.624 

80Z0 

88Z0 

6980 

0.925 

1.000 

aw% 

x" 


c 

c 

c 

c 

c 

c 

00 

C\ 

T- 

c 

cc 

c 

c\ 

c 

On 

CD 

CO 

C 

0 407 

069 0 

608 0 

000 1- 



%.4/i J) — average ahsohile dev i at ion *“ defined as n/).v((exp calc)! exp) 1 OO)/ hi 



Table 4.22: Vapor-iiquid equilibrium pressures and phase compositions for R-134a(1)R-141b(2} mixtures 
with different models. 

j 

P(MPa) 


UNIQUAC-R 

I 

0.245 

0.356 

0.468 

0.550 

0.652 

0.782 

0.894 

0.974 

1.064 

1.232 

1.379 

1.481 

1.680 

1.026 1 

o 

(0 

o 

UNIQUAC 


0.245 

0.356 

0.468 

0.550 

CO 

to 

CO 

d 

0.784 

9680 

ZZ60 

1.067 

1.234 

1.380 

1.482 

1.681 

1.016 


NRTL 


0.245 

0.354 

CD 

d 

0.546 

0.649 

08Z0 

0.894 

9Z6 0 

1.067 

1.236 

1.382 

1.483 

1.684 

0.874 


expti 


m 

, 

CO 

CO 

CO 

CD 

CM 

CD 

0.340 

0.461 

0990 

0.649 

98Z0 

CO 

CO 

00 

d 

9Z60 

CD 

CD 

O 

CD 

CO 

CM 

1.389 

1.504 

1.685 

%AAD 

[ 

calc 

UNIQUAC-R 

O 

O 

O 

d 

CM 

CO 

d 

L 

0.503 

6890 

Z990 

0.739 

98Z0 

0.815 

CO 

CO 

d 

CO 

CO 

d 

i 

0.925 

0.949 

1.000 

3.759 

UNIQUAC 

O 

O 

O 

d 

0,325 

£090 

6890 

0.667 

o 

d 

Z8Z0 

0.815 

CO 

00 

d 

6880 

0.925 

0960 

o 

o 

o 

CD 

CO 

CO 

NRTL 

o 

o 

o 

d 

0.320 

0090 

Z890 

Z990 

0.741 

0.789 

CD 

d 

CD 

CO 

d 

0.891 

0.926 

0960 

1.000 

3.731 

expti 


o 

o 

o 

d 

0.278 

0.479 

0.576 

0.643 

0.712 

69Z 0 

0.793 

0.822 

8Z8 0 

0.917 

0.944 

000 i. 

%AAD 

x" 


o 

c 

c 

c 

0.056 

n 117 

• 

• (£ 

i C 

n 226 

1 

1 T- 

l CC 

! c 

o 

a: 

CC 

c 

0 450 

> c 
i c\ 

‘ ID 

! c 

) N- 
i iC 
) CD 
) C 

0 777 

. N- 
. LD 
. CO 

1 c 

1 nob 

1 

i 

1 


%/l A D avercii^e absolute cieviat ion — defined as a/).y ((exj-) calc ) / exp ) * 1 OO)/ N 




rabfe 4.23: Vapor-liquid equilibrium pr^sures and phase compositions for R-23{1)/R-13(2) mixtures 
with different models (contd.). 

_ — — — 1 

P(MPa) 


UNIQUAC-R 


1.922 

2.090 

2.235 

2.390 

o 

V— 

LO 

CM 

5 

CM 

2.682 

2.708 

2./12 

2./14 

2.709 

O 

o 

04 

2.618 

2.554 

2.429 

2.309 

•1.549 1 

o 

03 

O 

NRTL UNIQUAC 


1.949 

2.119 

2.266 

2.423 

2.545 

2.650 

2.719 

2.746 

2.750 

2.751 

2.746 

2.712 

2.654 

2.589 

2.463 

2.341 

1.549 

. 

1.949 

2.103 

2.241 

2.389 

2.507 

2.619 

2.709 

, 2.753 

2.763 

2.768 

2.774 

2.751 

2.704 

2.645 

2.511 

2.341 

1.324 


expti 


272.94 K 

1.967 

2.135 

2.279 

2.423 

2.529 

2.625 

2.694 

2.733 

2.749 

2.752 

2.762 

2.751 

2.723 

2.683 

2.597 

2.499 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

0.107 

0.199 

CO 

C73 

CM 

d 

0.378 

0.455 

CM 

LO 

d 

0.564 

6Z90 

6890 

o 

CO 

CD 

d 

0690 

0.752 

6080 

6060 

1.000 

3.570 

UNIQUAC 

o 

o 

o 

d 

o 

d 

66 rO 

0.298 

CO 

CO 

d 

0.455 

0.521 

0.564 

6Z90 

6890 

o 

CO 

CO 

d 

o 

CD 

CD 

d 

0.752 

6080 

606 0 

O 

o 

o 

3.564 

NRTL 


o 

o 

o 

d 

o 

o 

d 

0.185 

0.274 

CO 

CO 

d 

0.425 

0.504 

0.562 

0.582 

968 0 

CO 

CD 

d 

0.713 

69Z0 

CO 

T“ 

CO 

d 

0.904 

o 

o 

q 

1.124 

expti 


o 

o 

o 

o 

860 0 

0.189 

0.278 

0.348 

0.436 

9080 

0.564 

0.586 

869 0 

0.652 

0.724 

0.783 

0.836 

0.919 

000 V 

aw% 

x“ 


c 

c 

c 

c 

c\ 

cc 

c 

c 

0 124 

c\ 
I a 
c\ 

c 

h- 

c\ 

c 

CC 

cc 

c 

c\ 

cc 

XT 

c 

oc 

CO 

LC 

c 

CC 

iC 

c 

c\ 

oc 

uo 

c 

n 652 

h- 

CO 

j c 

CC 

c 

cc 

c 

n 860 

CC 

cc 

a 

c 

C 

C 

c 



ahsohUe dcvialion - defined as «/7.v((exp — calc)! exj)) * 1 OO)/ N 


Table 4.23: VapoNiquid equilibrium pressures and phase compositions for R-23(1)/R-13(2) mixtures 
with different g'^ models (contd.). 



P(MPa) 

calc 

UNIQUAC I UNIQUAC-R 


1.005 

1.152 

1.285 

1.394 

1.479 

1.543 

1.619 

1.660 

1.691 

1.707 

1.714 

1.705 

1.680 

1.648 

o 

CD 

1.542 

CD 

CD 

CO 

CO 

4.554 


1.018 

00 

CD 

CO 

o 

CO 

CO 

1.500 

999 T 

CM 

CD 

CO 

CO 

CD 

1.715 

1.730 

1.738 

1.729 

1.704 

1.671 

1.624 

1.564 

1.486 

CO 

4.532 

NRTL 


1.018 

1.177 

CO 

CO 

1.420 

1.502 

1.563 

1.635 

h- 

CD 

1.704 

00 

1.724 

1.714 

1.688 

1.655 

1.607 

1.546 

1.467 

CO 

T— 

T-* 

4.420 

expti 


is:: 

CD 

CD 

CD 

CN 

! 

1.218 

1.306 

CD 

00 

CO 

1.458 

1.518 

1.566 

1.623 

1.654 

1.677 

1.692 

1.701 

1.703 

1.691 

1.672 

1.639 

1.606 

1.574 

1.499 

%AAD 


calc 

UNIQUAC-R 

0000 

CO 

CO 

CD 

0.231 

0.302 

0.356 

0.397 

0.451 

0.487 

0.523 

0.551 

0690 

0.623 

Z990 

0.684 

d 

CD 

d 

0.787 

6660 

12.740 

UNIQUAC 

O 

o 

o 

CD 

CO 

CO 

o 

0.231 

0.302 

0.356 

0.398 

0.452 

0.487 

0.523 

0.551 

0690 

0.623 

0.657 

n 

CO 

CD 

d 

00 

d 

0.746 

0.786 

6660 

12.817 

NRTL 

O 

O 

o 

CD 

CO 

CD 

CO 

CO 

CM 

d 

0.301 

0.352 

0.392 

0.445 

0.481 

00 

d 

0.547 

0.587 

0.622 

0.657 

9890 

CD 

d 

09Z0 

0.793 

6660 

12.699 

expti 


O 

O 

O 

CD 

CD 

CO 

O 

CD 

0.167 

0.237 

0.296 

CO 

d 

0.429 

0.475 

0.514 

0.549 

909 0 

9990 

0.712 

0.764 

0.821 

Z980 

606 0 

O 

O 

q 

aw% 

x" 


o 

o 

o 

o 

n 043 

0.092 

0.146 

0.200 

0253 

0.336 

0 397 

0.463 

0.516 

CD 

LC 

C 

0 663 

o: 

oc 

r-' 

> c 

008 0 

0 859 

n 903 

0 938 

c 

c 

c 



%AAD - average ahsohite deviation- defined ay (^a/>.v((exp-6'(:r/6')/exp)*10o)/ N 



Table 4.23: Vapor-liquid equilibrium pressures and phase compositions for R-23(1)/R-13(2) mixtures 
with different models (contd.). 

— — — 1 

P(MPa) 


UNIQUAC-R 

I 

0.435 

0.473 

0.509 

0.534 

0.566 

0.587 

0.604 

0.614 

0.617 

0.617 

0.612 

9090 

0.596 

0.577 

0.557 

0.530 

0.509 

0.474 

1.738 1 

calc 

UNIQUAC 


0.441 

0.480 

0.516 

0.542 

0.573 

9650 

0.612 

0.623 

0.625 

0.626 

0.620 

0.615 

0.604 

5850 

595 0 

0.537 

0.516 

0.480 

1.737 


NRTL 


0.427 

0.471 

0.510 

0.536 

0.568 

065 0 

Z09 0 

0.621 

0.624 

0.625 

2290 

0.617 

CO 

o 

CD 

d 

0.592 

0.575 

CO 

ID 

d 

0.525 

0.480 

1.603 

expti 


S 

Csi 

CN 

CD 

d 

CO 

CD 

d 

0.520 

0.543 

0Z90 

6850 

5090 

0.617 

o 

CM 

CD 

d 

0.622 

o 

CM 

CD 

d 

0.617 

o 

CD 

d 

Z650 

0.583 

0.562 

0.546 

0.513 

%AAD 


calc 

UNIQUAC-R 

o 

o 

o 

d 

0.109 

0.203 

0.266 

CO 

CO 

d 

0.399 

0.449 

0.495 

0.510 

0.521 

0.590 

0.613 

0.650 

0.707 

o 

CD 

d 

0.835 

0.894 

1.000 

O 

CO 

M- 

CO 

UNIQUAC 
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Table 4.23: Vapor-liquid equilibrium pressures and phase compositions for R-23(1)/R-13{2) mixtures 

with different models, , 
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g.4.1 5 Vapor-liquid equilibrium pressures and phase compositions for 
-23(1 )/R-1 3(2) mixtures at different temperatures (K). 
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4.24: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-123(2) mixtures 
with different models. 
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Table 4.25: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-114(2) mixtures 
with different models (contd.). 
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Table 4.25: Vapor-liquid equilibrium pressures and phase compositions for R-22{1)/R-114(2) mixtures 
with different models (contd.). 
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bla 4.25: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-114(2) mixtures 
with different g'^ models. 
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Tabie 4.26: Vapor-liquid equilibrium pressures and phase compositions for R-1 2(1 )/R-1 14(2) mixtures 
with different models. 
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Fig.4.18 Vapor-iiquid equilibrium pressures and phase compositions for 
R-1 2(1 )/R-1 14(2) mixtures at different temperatures (K). 





rabie 4.27: Vapor-liquid equilibrium pressures and phase compositions for R13B1(1)/R-12(2) mixtures 
with different models (contd.). 
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Table 4.28:Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-142b(2) mixtures 
with different models (contd.). 
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able 4.28;Vapor-iiquid equilibrium pressures and phase compositions for R-22{1)/R-142b(2) mixtures 
with different models (contd.). 
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able 4.28:Vapor“liquid equilibrium pressures and phase compositions for R-22(1)/R-142b(2) mixtures 
V / h different models (contd.). 
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Table 4 . 28 :Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-142b{2) mixtures 
with different models. 
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4.29: Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-142b{2) mixtures 
with different models (contd.). 
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Table 4.29: Vapor-liquid equilibrium pressures and phase compositions for R-32{1)/R-142b(2) mixtures 
with different models. 
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Table 4.30:Vapor-liquicl equilibrium pressures and phase compositions for R-32(1)/R-123(2) mixtures 
with different models (contd.). 
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Table 4 . 30 :Vapor-liquid equilibrium pressures and phase compositions for R-32(1)/R-123{2) mixtures 
with different models. 
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Table 4.31: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-600a(2) mixtures 
with different models (contd.). 
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Table 4.31: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-600a(2) mixtures 
with different models. 
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able 4.33: Vapor-liquid equilibrium pressures and phase compositions for R-600a{1)/R-236fa(2) mixtures at 303.68 K 
with different models. 


P(MPa) 

calc 

UNIQUAC 1 UNIQUAC-R 


0.325 

0.389 

0.437 

0.491 

ZOSO 

0.522 

0,530 

0.536 

0.537 

0.535 

0.531 

t — 

LD 

d 

0.476 

0.440 

0.410 

6080 ■ 


0,325 

6800 

CD 

CO 

d 

0.491 

0.507 

0.522 

0.530 

0.536 

0.537 

0.535 

0.530 

0.512 

0.476 

0.440 

0.410 

0.810 


NRTL 


0.325 

0.389 

00 

d 

CO 

CO 

d 

0,498 

0.514 

0.524 

0.534 

0.536 

0.535 

0.532 

0.517 

0.487 

0.449 

0.410 

0,237 


expti 



0.326 

0.390 

0.437 

0,483 

CO 

CD 

d 

0.514 

0.523 

CO 

CO 

ID 

d 

0.535 

0.535 

0.531 

0.517 

oo 

d 

0.449 

0,411 

%AAD 


calc 

UNIQUAC-R 

] 

o 

o 

o 

d 

0.199 

CD 

CO 

■ 

o 

0.437 

0.472 

0.508 

o 

CO 

ID 

d 

0.555 

999 0 

0.595 

0.617 

CO 

CD 

d 

0.791 

. 0060 

o 

o 

p 

3.206 

UNIQUAC 

o 

o 

o 

d 

0.198 

0.316 

0.437 

0.472 

8090 

0.530 

0.555 

999 0 

969 0 

0,617 

0.684 

0,791 

006 0 

o 

o 

o 

3.012 

NRTL 

0000 

0.194 

o 

CO 

d 

0.415 

0.449 

0.488 

0.515 

0.556 

0.578 

0.620 

0.643 

OOZ'O 

98Z0 

0.887 

1.000 

9360 

expti 


O 

O 

O 

d 

0.201 

0.302 

CO 

o 

d 

CO 

d 

CN 

CO 

d 

0.510 

0.553 

ZZ90 

0.624 

tD 

CD 

d 

o 

d 

0.791 

0.892 

o 

o 

o 

%AAD 



o 

o 

o 

o 

0.061 

0.120 

0.219 

0.265 

0.330 

0.388 

0.490 

0 551 

£99 0 

0.7l'5 

0.810 

868 0 

896 0 

1.000 



o/i/l / ) = avcr(i}^e absolute deviation — defined as ahs{{exp calc) / exp ) * 1 oo)/ N 










Table 4.34: Vapor-liquid equilibrium pressures and phase compositions for R-22{1)/R-23{2) mixtures 
with different models (contd.). 


P(MPa) 

calc 

UNIQUAC-R 


CO 

00 

o 

CM 

O 

h- 

d 

o 

CD 

d 

1.240 

1.556 

1.799 

2.003 


N- 

CO 

CM 

CM 

00 

CO 

MT 

CM 

o 

CD 

CO 

CM 


9680 

1.139 

CO 

o 

CD 

— 


LO 

CD 

M" 

CM 

3.006 

3.749 

3.817 

4.115 

1.720 

UNIQUAC 

m 

CO 

N. 

CN 

0.488 

0.702 

to 

o 

CD 

d 

1.240 

CO 

to 

1.799 

2.003 

2.126 

D- 

CO 

CM 

CM 

uo 

CO 

M- 

CM 

O 

CD 

CO 

CM 



1.108 

o 

cq 

t — 

2.000 

CD 

CM 

2.995 

3.745 

1 . 

3.814 

LO 

T** 

t — 

1.719 

NRTL 

0.489 



1 

LO 

S 

T— 

CO 

N- 

1.978 


CD 

CM 

CM 

CM 


r-- 

M- 

CM 

CM 



00 

O 

o 

LO 

2.000 

r^ 

CD 

CM 



3.814 

4.142 

1.192 

expti 


0.495 

0.715 

m 

o 

CD 

d 

uo 

o 

CM 

•t— 

LO 

UD 

1.800 

UO 

CO 

o 

CM 

2.205 

2.355 

O 

ID 

CM 

Q 

< 

< 


906 0 

o 

T— 

LO 

o 

lO 

2.000 

LO 

CD 

CM 

2.995 

3.500 

3.850 

4.123 

Q 

< 

< 

>; 

calc 

UNIQUAC-R 

o 

o 

o 

d 

0.346 

0.514 

ZZ90 

08Z0 

0.853 

906 0 

0.940 

0.969 

1.000 


293.15 K 

o 

o 

o 

d 

0.207 

0.455 

0.632 

0.744 

0.833 

0.841 

9960 

1,000 

3.375 

UNIQUAC 

o 

o 

o 

d 

T— 

CO 

d 

CD 

CM 

UD 

d 

CO 

o 

d 

D- 

O 

00 

d 

CD 

CD 

00 

d 

CM 

5 

d 

0.941 

00 

CO 

CD 

d 

1.000 

1.810 

o 

o 

o 

d 

CO 

T— 

CM 

d 

CO 

XJ- 

d 

00 

CM 

CO 

d 

lO 

LO 

d 


0.949 

8960 

1 .000 

3.375 

NRTL 

o 

o 

o 

d 

CO 

CD 

CO 

d 

05 

CM 

to 

d 




CO 

o 

CD 

d 

0.938 

CO 

CO 

CD 

d 

1.000 

1.129 

o 

o 

o 

d 

CO 

o 

CM 

d 

CM 

LO 

d 

o 

CO 

CO 

d 


00 

CO 

d 

0.947 

0.957 

1,000 

2.375 

expti 

o 

o 

o 

d 

0.346 



o 

CO 

d 

CO 

LO 

CO 

d 

CO 

o 

d 

0.940 

CD 

CD 

CD 

d 

o 

O 

O 

-T — 

%AAD 

o 

o 

o 

d 

0.207 





0.910 

9960 

O 

o 

o 

%AAD 

x“ 

0000 

0.102 










o 

o 

o 

d 

CD 

O 

d 

CD 

CD 

d 

CM 

CD 

CO 

d 

T~ 

CO 

in 

d 

0.693 

0.841 

0.928 

1.000 



D = average absolute deviation - defined as fl/?i ((exp - calc) / exp) * 1 
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o 

V. 

3 

X 


CN 

co" 

CM 


CM 

CM 
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VC 

c 

o 

’Eo 

o 

a 

E 

o 

o 

<D 

t/i 

03 

£ 

a 

"O 

c 

03 

tfi 

0 

\m 

3 

(/) 

"\ 

a ^ 

E i 

3 a> 
■= 

1 E 

g-- 

0 

■D c 

1- ^ 

CT o 

"7 

o ‘-5 

Q. il 
«5 .ti 
> 5 

CO 

o 

Xi 

fS 


CD 


Cl 

X 

0 



1.164 

1.575 

2.015 

2.500 

2.989 

3.564 

4.051 

4.552 

4.690 

1.214 


1.965 

2.443 

3.021 



3.562 

3.995 

4,587 

5.014 

5.194 

5.256 

0.235 


1.154 

1.575 

2.015 

2.500 

2.990 

3.564 

4.055 

4.568 

4.782 

1.021 


1.965 

2.419 

eooe 

3.560 

4.001 

— — 

4.597 

r^ 

T— 

q 

id 

5.187 

5.244 

CO 

CM 

CO 

d 


CD 

1.500 

O 

cn 

CD 

2.503 

2.997 

3.566 

4.053 

4.578 

4.788 

CO 

d 


1.943 

2.419 

CO 

o 

o 

cd 

3.560 

o 

o 

4.597 

5.017 

5.187 

5.244 

0.323 


00 

00 

1.525 

2.020 

2.510 

2.990 

3.515 

o 

o 

o 

4.500 

4.805 

%AAD 


1.950 

2.430 

3.020 

3.560 

3.995 

4.575 

4.990 

5.170 

5.255 

%AAD 

303.15 K 

o 

o 

o 

d 

0.246 

0.467 

o 

(D 

d 

OOZO 

0.781 

0.851 

5060 

0.917 

1.425 

323.15 K 

O 

O 

o 

d 

0.215 

0.382 

0.491 

695 0 

0.627 

CD 

CD 

d 

(D 

CD 

d 

o 

00 

CD 

d 

2.693 

o 

o 

o 

d 

0.250 

0.449 

0.598 

0.702 

0.792 

0.852 

5060 

0.917 

1.563 

0000 

0.215 

0.382 

d 

6550 

0.630 

Z990 

SZ90 

o 

CD 

q 

o 

2.629 

000 0 

0.236 

CD 

d 

0.612 

o 

d 

0.794 

0.852 

0.911 

0.919 

2.046 

O 

O 

o 

d 

0.212 

0.384 

0.495 

0.561 

0.630 

5990 

0.677 

9Z90 

2.649 

O 

O 

O 

d 

0.246 

0.467 

to 

0 

ld 

1 

OOZO 

0.781 

0.845 

0.904 

0.936 

%AAD 


o 

o 

o 

d 

0 222 

0.396 

0.505 

0.577 

CD 

O 

0.682 

0.693 

0.691 

avv% 

o 

o 

o 

o 

060 0 

0.228 

0.367 

0496 

0.639 

9910 

s 

CO 

o 

0 922 



c 

c 

c 

c 

CC 

CT 

C 

c 

0 218 

N- 

Cn 

CC 

c 

CN 

c 

C\ 

IT 

C 

0 604 

0.641 

099 0 



%//,4 /-> = a\’cra}^a ahsohile ckvkUion — liafined as a/ys{{eKp — ca/a) / exp) I oo)/ 








Table 4.35:Vapor-liquid equilibrium pressures and phase compositions for R-152a(1)/R-150a(2) mixtures 
with different models (contd.). 


P(MPa) 


UNIQUAC-R 


0.037 

0.135 

0.211 

0.278 

0.347 

0.445 

0.490 

0.543 

0.594 

0.642 

989 0 

0.722 

1.990 1 

o 

cn 

o 

UNIQUAC 


0.037 

0.134 

0.207 

0.271 

0.336 

CO 

d 

0.477 

0.532 

0.587 

0.638 

9890 

0.721 

88Z0 


NRTL 


0.037 

CD 

CO 

d 

CO 

o 

CM 

d 

0.271 

CD 

CO 

CO 

d 

0.431 

0.477 

0.533 

CO 

CO 

CO 

d 

0.639 

CO 

00 

CO 

d 

0.722 

1.054 

exptl 


<N 

CO 

O 

CO 

zeoo 

CO 

d 

0.205 

0.281 

CD 

CO 

CO 

d 

0.431 

0.470 

0.530 

0.595 

0.639 

9890 

0.721 

%AAD 


calc 

UNIQUAC-R 

O 

o 

o 

d 

0.743 

0.846 

0.891 

0.920 

0.948 

8960 

0.969 

6Z60 

Z860 

0.994 

1.000 

2.398 

UNIQUAC 

o 

o 

o 

d 

o 

d 

\ 

0.842 

CO 

CO 

CO 

d 

0.917 

0.945 

996 0 

Z960 

ZZ60 

CD 

00 

<J) 

d 

0.994 

o 

o 

o 

2.614 

NRTL 

o 

o 

o 

d 

0.745 

0.844 

00 

00 

00 

d 

0.918 

0.946 

0.957 

CO 

CD 

CD 

d 

CO 

(j) 

d 

Z860 

0.994 

o 

o 

o 

2.776 

exptl 


o 

o 

o 

d 

0.693 

0.803 

o 

CD 

CO 

d 

0.891 

0.927 

0.940 

8960 

996 0 

0.977 

986 0 

1.000 

%AAD 

x" 


o 

o 

o 

d 

0 106 

0.194 

0 278 

0.372 

cr: 

Cs 

LC 

c 

809 0 

0.705 

0799 ^ 

0.881 

c 

ic 

a 

c 

000 1 








o 

3 

X 

S 

'fC 
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to 

T“ 

(k 


Ct3 

04 

lO 

T“ 

I 

Of 

L« 

a. 

m 

c 

o 

"35 

o 
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E 

o 

o 

o 

0) 

cc 

x: 

a 

X5 

c 

CO 

0) 

o 

3 

C/) 

(/) 

<D 


E i2 

3 O 
*2 

1 1 

3 

Q cn 
"5 c 

3 £ 

cr o 
’•T !J= 
O ’"5 
a x: 
ns 

> 5 

li) 

CO 

Q 

XS 

CO 

H 


cr 

1 

o 

< 

3 

a 


CL 

X 

0 



OBOO 

0.202 

0.281 

0.417 

0.503 

0.681 

d 

0.825 

0.971 

1.036 

1.152 

1.229 

4.424 1 


6Z00 

0.201 

0.277 

00 

o 

d 

o 

G 

d 

6S90 

80Z0 

008 0 

096 0 

1,020 

1.149 

1.230 

4.353 


6Z00 

0.214 

0.293 

0.423 

0.501 

0.661 

80Z0 

0.797 

0.947 

1.019 

1,149 

1.230 

5.115 

/ 

4 

i 

4 

O 

00 

o 

o 

0.158 

0.277 

0.440 

0.515 

0Z90 

80Z0 

0.791 

0.942 

1.016 

■r— 

1.229 

%AAD 

o 

o 

o 

d 

0.637 

0.753 

00 

d 

CD 

00 

00 

d 

CO 

G 

d 

0.940 

SS6 0 

0.973 

0860 

0.992 

1.000 

m 

CD 

o 

o 

o 

d 

0.634 

h- 

d 

r^ 

OO 

d 

\ 

3880 

0.927 

0.937 

0.952 

0.971 

0.979 

0.992 

1.000 

T— 

00 

G 

o 

o 

o 

d 

0.652 

0.757 

0.845 

ZZ80 

0.921 

0.930 

0.946 

Z96 0 

0.976 

0.991 

1 ,000 

1.786 

o 

o 

o 

d 

0.634 

89Z0 

0.825 

0.855 

0.902 

't — 

G 

d 

0.928 

0.952 

0.964 

0.982 

1.000 

%AAD 

o 

c 

c 

G 

0102 

0.164 

0 272 

OC 

c 

0.488 

c\ 

cc 

ic 

c 

CC 

T- 

<£ 

C 

T- 

G 

c 

0 827 

G 

CC 

a 

c 

C 
• C 

c 

: T- 








absolute deviation — defined as a/?.v((exj3 cv/Zc*) / exp) * 1 OO)/ A/ 






Table 4.36: Vapor-liquid equilibrium pressures and phase compositions for R-134a{1)/R-142b(2) mixtures 
with different models (contd.). 


0^ 

O) 

CN 


N- 

o 

00 

CD 

N. 

CN 

CN 

T-" 


CD 


00 

o 

T— 

CN 

CO 
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T— 



T— 

T— 

CN 

CN 

CN 

CN 

; CN 

CN 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 


jj 
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CD 

CD 

CN 

UO 

N. 

CD 

CO 

CD 

r^ 

xr— 

CO 

cu 

o 

o 




CD 


00 

CD 

xr- 

CN 

CO 

'N* 

CD 

2 
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x— 

T— 

CN 

CN 

CN 
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d 

d 

d 

d 

d 

d 
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m 

0 

u, 

3 

X 

1 

21 

n 

CM 


CO 

CO 

T“ 

I 

a: 


</) 
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o 


CO 

o 

a 
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o 

o 

o 

CO 

CO 

sz 

CL 

“D 
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CO 

CO 

o 

3 

<0 

CO 
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a 

E i 

3 O 

*5 

Xi O 

E £ 

g-u. 

o 


O) 


CD 

CO 

,Q 

m 

H 


P(MPa) 


UNIQUAC-R 


0.337 

0.3/2 

0.404 

d 

0.4/2 

O.bOb 

0.539 

0.578 

n 

ID 

d 

0.643 

0990 

0.141 1 

o 

CO 

o 

UNIQUAC 


0.337 

0.372 

0.404 

0.439 

0.472 

9090 

0.539 

8Z90 

0.615 

0.643 

0990 

O 

d 


NRTL 


0.337 

CO 

CO 

d 

0.406 

0.441 

CO 

d 

0.505 

0.538 
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s 

d 

CO 

CD 

d 
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608 0 
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o 

00 

cn 
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CO 

CO 

CD 
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>; 
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CM 

CD 
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4.37: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-12(2) mixtures 
with different models. 
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Table 4.38: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-152a{2) mixtures 
with different models. 
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Table 4.40: Vapor-liquid equilibrium pressures and phase compositions for R-134a(1)/R-152a(2) mixtures at 199.80 K 
with different models. 
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Table 4.41: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-134a(2) mixtures 
with different models (contd.). 
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fable 4.41: Vapor-liquid equilibrium pressures and phase compositions for R-22(1)/R-134a(2) mixtures 
with different models. 
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through 4.7 provide comparisons of the fitted WS model using the NRTL equation with 
the measured total pressure for these systems. 

Tables 4.9 through 4.41 give the VLB results for 32 systems for which P-x-y data were 
available; calculated pressures and vapor mole fractions are compared with the respective 
experimental values and overall average absolute deviations (%AAD) are presented. 
Figures 4.8 through 4.27 provide comparisons of the fitted WS model using the NRTL 
equation with the measured total pressure and vapor composition for some selected 
systems included in Table 4.9 through 4.41. 

The present study reveals that the WS model using any of the three equations allows, 
in general, accurate representation of phase behavior for a variety of refrigerant mixtures 
over a range of temperatures and pressures . Moreover, it is interesting to note that the 
results with the UNIQUAC-R model are comparable to those obtained with the 
UNIQUAC model, signifying that the combinatorial term has little influence on the 
accuracy of VLB predictions over the pressure range of interest. 

4.4 Conclusions 

The results showed that the Wong-Sandler mixing rules with the NRTL, UNIQUAC or 
UNIQUAC-R activity coefficient models in all cases could be accurately fitted to binary 
refrigerant VLB data. However, when using thislnodel it is important that the appropriate 
initial guesses be made for the mixing rule parameter, ki 2 , and the model parameters, 
'^ 21 - Although the parameters in the models are slightly temperature dependent, 
reasonably good predictions can be achieved at other temperatures even if these 
parameters are assumed to be temperature independent. This would allow extrapolation to 
higher and lower temperatures when a very limited amount of experimental data are 
available for a system. Where VLB data exist for more than one temperature, 
simultaneously fitting the data at all temperatures should provide the best parameters to 
use over the entire temperature range. Additional calculations and predictions should be 



made to further verify this conclusion. Another alternative would be to simultaneously fit 
VLB along one isotherm and excess enthalpy data, if such data were available. 

In summary, on the basis of the results presented here, the WS model appears to be 
promising in correlation and extrapolation of the vapor-liquid equilibria of refidgerant 
mixtures over a range of temperatures and pressures, especially for the newer and 
unconventional mixtures now imder consideration. 
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